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Abstract 

We study a refined large spin limit for twist operators in the sl(2) sector of AdS/CFT. 
We derive a novel non-perturbative equation for the generalized two-parameter scaling 
function associated to this limit, and analyze it at weak coupling. It is expected to smoothly 
interpolate between weakly coupled gauge theory and string theory at strong coupling. 



1 Introduction, Main Result, and Open Problem 



The perhaps most interesting subset of all local composite quantum operators of planar 
M = 4 supersymmetric gauge theory is formed by the sector of sl(2) twist operators. The 
reason is that these bear many similarities with the twist operators of QCD. They may be 
symbolically written as 

Tr (V M Z L ) + . . . , (1.1) 

which is a shorthand notation for intricate linear superpositions of all states where the M 
covariant derivatives D act in all possible ways on the L complex scalar fields Z. Here L is 
a su(4) R-charge, frequently denoted as J in the literature, and M is a Lorentz spin, often 
called S. Our labelling refers to the magnetic spin chain picture of these operators, where 
L is the length of the chain, and M is the "magnon number" . The twist of an operator is 
defined as the classical dimension minus its Lorentz spin, so the length L equals the twist 
in the case of (II. ip . 

M = 4 gauge theory is a superconformal field theory. Therefore proper superpositions 
of the operators (11.11) must carry a definite charge A under dilatations. It generically 
is, in contradistinction to the R-charge L and the Lorentz charge M, coupling constant 
dependent: A = A(g). Its anomalous part j(g) is defined ad3 

A(g)=M + L + 1 (g), (1.2) 

where M + L is the classical dimension of the operators (II. ip . In the case of the operators 
(II. ip 'j(g) behaves in a very interesting way as the spin M gets large at fixed twist L. It 
grows logarithmically with M at all orders of the coupling constant g defined as 

9 > = §^ * (1.3) 

7T Z 16 7T 

where A is the 't Hooft coupling. The prefactor of the logarithm is a function of g. We call 
it the universal scaling function f(g): 

A - M - L = 7 G?) = f{g) logM + . . . . (1.4) 

This behavior is a special case of so-called Sudakov scaling, see (TJ . In the twist L = 2 
case it equals twice the cusp anomalous dimension of light-like Wilson loops [2J. The 
independence or "universality" of the function f(g) on the twist L, with L arbitrary but 
finite, or even L — > oo as long as L «C logM, was first pointed out at one loop in [3j, and 
conjectured to hold at arbitrary loop order in [6]. It would be very interesting to rigorously 
prove that the twist-two (L = 2, M — > oo), twist- L (L fixed, M — > oo), and the universal 
(L, M — > oo,L logM) scaling functions f(g) of the operators (11.11) indeed all coincide 
for arbitrary values of g: f(g) = fW(g) = f^(g) = ■■■ = f^ L \g) = ... = f^\g). 



1 The anomalous dimension j(g) is related to the energy E{g) of the integrable long range spin chain 
describing the operators through j(g) = 2g 2 E(g). It should not be confused with the energy of 

string states which equals A(g) via the AdS/CFT correspondence. 
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Based on the conjectured all-loop integrability of planar M = 4 theory [I], the weak 
coupling expansion of f(g) is known from the solution of an integral equation obtained 
from the asymptotic Bethe ansatz for these operators El [7] . It agrees to four order^l 
with field theory [8]: 

/(<?) = 8 g 2 - f TrY + | A 6 - 16 (J^ ^ + 4 C(3) 2 ) g* ± . . . . (1.5) 

Testing the Bethe ansatz to five orders in field theory might not be out of reach [9]. 

The strong coupling expansion may also be obtained from the same integral equation [7j 
which generates the small g expansion (jl.5p . After the initial studies PHtH], an impressive 
analytical expansion method to any desired order was worked out in [T2j. The starting 
point of this systematic approach was an important decoupling method discovered by Eden 
[13] . see also pH]. The series starts as 

3 log 2 K 1 

f(g)=4g — . . . , 1.6 

ii ATI 1 g 

where K= j3(2) = Catalan's constant. The first two terms on the r.h.s. agree, respectively, 
with the classical and one-loop [T4"l [15] result from semi-classical string theory, and the 
third term is the two- loop correction very recently obtained in [U, [18]. It would be very 
interesting to also check the three-loop term in string theorjIE 

So it appears that f(g) is the first example of an exactly known, via the solution 
of a linear integral equation [7J, function which smoothly interpolates between a gauge 
theory and a string theory observable in the AdS/CFT system. A natural question is 
whether further interesting examples may be found, and whether the function f(g) may 
be generalized. A major obstacle is the fact that we currently only know the asymptotic 
spectrum of the planar M = 4 model, as was recently unequivocally established in [T9] . 
Important clues come from both taking a closer look at the scaling law (11.41) in the one-loop 
gauge theory [3] , and at intriguing string theory results [2H 121] generalizing the expansion 
(11.61) . Put together, these suggest that at weak coupling an interesting generalized scaling 
limit might exist, where3 

M — > oo , L — > oo , with j := — = fixed . (1.7) 

log M 

We will prove in this paper that this is indeed the case, first at one-loop order, and then 
beyond. More precisely, we will show that a generalized scaling function f(g,j) exists to 



2 To be precise, the four-loop field-theory result of [5] agrees numerically with the analytic Bethe ansatz 
prediction in (jl.5|) to 0.001%. An analytic proof would be most welcome. 

3 However, one important fact to keep in mind is that the weak coupling expansion (|1.5[) has a finite 
radius of convergence, while the strong coupling series (|1.6j) is asymptotic, and, apparently, not even Borel- 
summable |12j . So (|1.6p follows from knowing all terms in (jl.5| . but, conversely, knowing all terms of the 
string expansion (|1.6|l does not allow to reconstruct the gauge-theoretic perturbation series (jl.5| without 
further input. Unfortunately, it is currently not even known what the nature of this input might be. 

4 The variable j was first explicitly introduced (up to a factor of 1/2) in eq. (3.1) of [21 . 
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all orders in perturbation theory 

A — M — L = j(g) = f(g, j) log M + . . . , (1.8) 

where f(g,0) = f(g) in (11.41) . This extends the one-loop results in [3], and the all-loop 
result at j = of [HI Ej- The latter is possible since in the limit (11.71) L — > oo, and we may 
therefore use the asymptotic Bethe ansatz methodology of [5j El [7j. 

The final result of our analysis, presented in detail in the ensuing chapters, is the 
following integral equation 

a(t) = -J—^ (t{t, 0) - 4 jf dt' JC(t,t') a(t')j . (1.9) 

It is essentially identical in form to the "ES" (no dressing phase) [0] and "BES" (with 
proper dressing phase) [7j equations. The kernel corresponding to the generalized scaling 
limit is quite involved as it contains various contributions. It reads 

,?./ #\ 9 - / ,x T -, / , , Jo(2qt) sin at' t' 

K{t, = 2 2 ^(2^, 2^') + tf h t, t'; a) - -^-^ — — 

t Z 7T t 

poo 

-Ag 2 / dt"t" k{2gt,2gt")K h {t",t';a). (1.10) 

Here 

£(t, t') = £ (t, + + #d(*, (l.H) 
is the kernel of the "BES" equation, where 

f> ( ,,,, t.h{t)Ut')-t' Jq (t) J! (tQ 

Ko{t,t) = ^— ^ , (1.12) 

t'Ji(t) Mt')-tJ (t) .hit') 

Kl{t,t) = 12 _ t , 2 > ( L13 ) 

and the kernel encoding the effects of the dressing phase is given by the convolution 

poo ill 

K d {t,t') = 8g 2 dfR^gt'^^—k^gt"^') . (1.14) 
Jo e — i 

For a possible mechanism generating this type of convolution structure see [22]. The novel 
contributions generated by a non-vanishing j are encoded in the kernel 



, 1 t t cos( at') sm( at) — t' cos( at) sm( at') 
Kh(t, f; a) = ^ e-s ^ 1 t ,_ f2 1 ^ — e- , (1.15) 



as well as the explicit, rightmost term of the first line of ( 11.101) . and the further convolution 
in the second line of that equation. The index h of K^(t, t'; a) stands for "hole" , its meaning 
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will become clear below. The corrections of the refined limit depend on a "gap" parameter 
a whose interpretation will also be explained. Its relation to j is fixed by the constraint 

4 a 16 f°° , , . t sin at 

3 = / dta(t)e*^. (1.16) 

f Wo t 

Lastly, the generalized scaling function of fll.Sp is given by 

f(g,j)=j + 16a(0). (1.17) 

It is determined by first solving the integral equation (jl.9p with the kernel (ll.lOp for 
the fluctuation density a(t) = a(t; g, a) as a function of g and a. Then a is found as a 
function of j by inverting the relation (11.161) . i.e. by computing a(j). This then yields 
a(t) = ait; g, a(j)) as a function of g and j, and the generalized scaling function f(g,j) is 
finally obtained by evaluating the latter at t = 0, see ( 11. 17ft . 

As in [51 E] , in practise it appears impossible to produce a closed-form solution of the 
equation (II. 9p . In fact, we did not even find an explicit solution at one-loop order, i.e. for 
g = 0. It is however possible to solve it iteratively in a double-expansion in small g and 
small j. Excitingly, the obtained function appears to be "bi-analytic", i.e. analytic in g 
around g — at arbitrary finite values of j, and vice versa. We therefore believe that 
our equations actually hold for arbitrary values of g and j. The beginning of this double 
expansion may be found in (14.1 71) . (T4~T9l . (14 . 2 (J|) . (14 . 2 1 p . which we have displayed by giving 
the four-loop result of the functions fi(g) ■ ■ ■ f^g) defined through 



f(g,j) = f(g) + J2f^)j n - (1-18) 



n=l 



Our truncation at four loops 0(g 8 ) and C(j 4 ) is due to space limitations, and one easily 
generates many more orders in g 2 and j if needed. A curious fact is the absence of any 
terms of 0(j 2 ), i.e. the function /2(fi0 is zero. We will come back to this point shortly. 

A very interesting question is how f(g,j) behaves at strong coupling. Indeed we would 
like to make contact with the already known results from string theory [T5 | \16 \ l3| [20 ], |2T | 123] . 
A potential trouble is that in the semi-classical computations pioneered by Frolov and 
Tseytlin [TS] the coupling constant g is intricately entangled with the, respectively, AdS 5 
and S 5 charges M and L. In |J] the strong coupling limit of the dimension A of the 
operators ( 11. ip was predicted from the results of [151 [IB] on the energy of a folded string 
soliton (see also the discussion in [20l [2TJ [23], E3] ) • The prediction reads 



A dassical = M + LVl + z 2 + ... , (1.19) 

where^ 

M^oo,L^oo,sf^oo, with M>L and fixing — , - , z := Ag l0g A , (1.20) 

9 9 L 



D The contemporaneously appearing work [23] uses the notation £ = 1/z and A = Ang. 
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and A is some sea le8 The result (QUI) was derived in [21] from the asymptotic Bethe 
ansatz [""J [7] with approximate strong coupling (AFS [2"5] only) dressing phase. While this 
constitutes an important check of the Bethe ansatz method, this had to work in the sense 
that the dressing phase [25J was extracted from the integrable structure of classical string 
theory [261 EZ] • What is missing is a derivation from a solution of the exact ansatz j5j [7] 
which interpolates between weak and strong coupling. Now it is tempting to identify, in 
view of (TOT) . 

4<7 , 
z = —. 1.21 

J 

Certainly the condition M 3> L with M, L — ► oo is satisfied in the weak-coupling limiting 
procedure (II. 7p . The more questionable assumptions in semi-classical string theory, as 
far as concerns extrapolating weak coupling results, are the fixation of M/g, L/g (see also 
footnotes on the previous page). Proceeding under this caveat we could then rewrite (11.191) 
as i 

■ \ 2 

J 



A (l _ 1(al -.U-L: (4<7yi+(^-) j ) \ot-M + ... . (1.22) 

If we now also expand in small j we find 

A classica i -M-L= \ 4g-j + ^- + G(j 4 ) | logM + ... . (1.23) 

The leading term 4g agrees with the first term in ( 11.61) . The one-loop string correction to 
(II . 191) was computed in 



Ai-ioop = -f= ; \ zVl + z 2 - (1 + 2z 2 ) log z + VI 



z 



Aoop ~V\VTT 

-z 2 + 2(1 + z 2 ) log(l + z 2 ) - (1 + 2z 2 ) log Vl + 2z 2 } . (1.24) 

Taking z — ► exo it produces the second term on the r.h.s. of ( 11.61) . If we were to again expand 
in small j via (11.211) we would find j 2 logj terms. The result f |1.24j) was fully derived in 
[21] from the asymptotic Bethe ansatz [7J E] with approximate strong coupling (AFS + HL 
[2*5] ) dressing phase. Once again, this is an important cross-check on the consistency of 
the extraction of the one-loop correction of the dressing phase from one-loop string theory 
[251 [29], see a l so the very recent derivation [30j, but does not answer the question how the 
dimensions of the gauge theory states in ( ll.ip "flow" to the energies of string theory states 
as the coupling increases. 



6 The scale A actually being used in the string theory calculations in [T~~J [16l OH [2Ql EH [23] seems to be 
somewhat unclear. Is the proper scale (1) A = 4 irg or (2) A = L or (3) A = 1? Since these calculations 
start from fixing M/g and L/g it would seem that they require either (1) or (2). At weak coupling we 
definitely have (3), as we are proving to all orders in this paper. Understanding the crossover of scales as 
one moves from weak to strong coupling, or vice versa, should be very interesting. 
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An interesting insight into the structure of further quantum corrections, i.e. two-loop 
and higher, to (11.191) . 01. 241) was obtained in in the limit z — > oo. In a paper contem- 
poraneous with ours [23], an impressive direct two-loop string calculation, in this limit, is 
performed which agrees with the results of [21]. However, Roiban and Tseytlin argue in 
[23] that after resumming infinitely many terms of the form j 2 \og k j all terms of the form 
j 2 might vanish. They furthermore noticed that some initial support for these considera- 
tions is provided by a fascinating and curious byproduct of our derivation: The function 
/2G7) in the expansion (11.181) is exactly zero! This suggests that extrapolation between the 
result at small g and the result at large g might indeed work out. 

Therefore an exciting open problem not addressed in this paper is to now solve our 
equations at strong coupling g — > 00 in order to see whether any of the above string results 
are reproduced, and whether the extrapolation works out. In fact, our derivation does not 
assume j to be small, so we are hopeful that under the identification (ll.2ip the full strong 
coupling expansion of f(g,j), i.e. (11. 191) . (11.241) and all further corrections, in generalization 
of the beautiful expansion of [T2] at j = 0, will be obtained. As already mentioned this is 
however not assured, as we might run into an order-of-limits problem, namely (11.71) versus 
( ll.20p . It would also be important to gain an understanding how the states corresponding 
to the generalized scaling function fit into the general classification of classical integrable 
curves [271 El] and their quantum fluctuations. It should be very interesting to see how 
the parameter j in ( II .7ft relates to the "filling fractions" of the classical curve. 

This paper is organized as follows. In section [2] we extend the study in [3] and take 
a close look at the fine-structure of the large-spin M anomalous dimensions of (II. ip at 
one-loop order. We derive our results both using traditional techniques as well as more 
sophisticated ones involving so-called non-linear integral (or also "Destri-DeVega" ) equa- 
tions, see [32] and references therein. In section [3] we generalize the methodology of the 
non-linear integral equations to all orders in the coupling constant and compute some novel 
finite size O(M ) corrections to the scaling behavior (11.41) . In section H] we extend our one- 
loop results to all loops, prove the existence of the novel generalized scaling function in 
( II. 8p . and derive the above equations determining it. 



2 One-Loop Theory 
2.1 Magnons and Holes 

The one-loop diagonalization problem of the operators (II. ip is equivalent to the one of an 
integrable spin chain with sl(2) symmetry. This was first discovered in [3"3"1 13"4"] and more 
specifically in the J\f — 4 context, extending the discovery of [35], in [36]. The allows to 
apply the Bethe ansatz, which then leads to the following one-loop Bethe equations 

\ a k 2 / j=l k j ' 
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where L is the length (=twist in this case) and M is the number of magnons, see (11.11) . 
The cyclicity constraint and the one- loop anomalous dimension 71 (see ( II. 2p ) are 



With the help of the Baxter function 

M 



M 



n^| = l and - n ^M =25:-^. (2.2) 



9=0 fc= l M fc t- 4 



Q(it) = " ( 2 - 3 ) 
fc=i 

one can write down an off-shell version of these equations 

m + Q(u + i)+(u-^\ Q(u-i) = t(u)Q(u), (2.4) 
where 

L 

t{u) =2u L + ^q i u L ~ i (2.5) 

i=2 

is the transfer matrix given in terms of the charges. The ground state for arbitrary L 
and M is unique and thus the corresponding charges are fixed. Clearly setting u = Uk 
in equation (12 .4p brings us back to (12. ip . However one of the advantages of (12 .4p is the 
possibility of identification of complementary solutions u = to (12.11) [3j. They are 
found as the zeros of the transfer matrix, i.e. from t(u) = and describe "holes". We thus 
have 

L 

t(u) = 2H(u-u ( £ ) ). (2.6) 

fc=i 

We can intuitively think of the hole roots as rapidities describing the motion of the Z- 
particles in the spin chain interpretation of the operators ( 11.11) . For a general value of L 

the equation (12. 6p has L solutions and thus there are L holes. One can prove that for any 

(k) 

state all magnon roots Uk and all hole roots are real. It is possible to find the q2 charge 
analytically by matching the three highest powers of u in the Baxter equation (12.41) : 

q 2 = --L(L—1)—LM — M (M — 1) . (2.7) 

Because q 2 and all higher charges explicitly depend on M the roots of t(u) will also, 
generically, depend on M. One can argue, however, that for the groundstate, and in the 
case L <C M, two of them are special, see Their magnitude is larger than the one of 
any other (hole or magnon) Bethe root, and scales with M as the magnon number M gets 
large, see [3]. To identify these roots one recalls that the mode numbers for magnons for 
the ground states, when L M, are given by [6] 

= k + — 7p— sgn(k) for k = ±1 ± 2, ±^ . (2.8) 
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The absolute value of the roots grows monotonically with \rik\- It follows from (I2.8P that 
the rapidities of the magnons and holes are parity- invariant. Among the holes there are 
two 'universal holes' which occupy the highest allowed mode numbers 

^L + M-l < = _L±M-1_ (2 . 9) 

The corresponding hole roots are precisely the one that scale with M. The remaining holes 
fill the gap in the mode numbers of magnons 

f L-3 L-3) 

<e • (2.10) 



2 2 

For the ground state, when L <C M, the magnitudes of the roots are thus ordered as 

\v£' 2) \>\u k \>v% ) (jV1,2). (2.11) 

2.2 The Counting Function and the NLIE 

A nice way to exploit the existence of the hidden hole degrees of freedom employs the 
so-called counting function, see [32J and references therein. It is defined as 

Z{u) = L<l>{u,\)+Y<t>{u-u k ,l) where <j){u, £) = i log ( V (2.12) 

Its name stems from the fact that, as one immediately sees from the definition (12.121) . 
Z(±oo) = ±7r(L + M) while the Bethe equations for the magnons and holes may be, 
respectively, expressed as 

Z(u j )=n(2n j + 5- 1) j = l,...,M, (2.13) 

Z(u^) = vr (2n[ k) + S - 1) fc = 1, . . . , L , (2.14) 

where 

S = L + M mod 2. (2.15) 

So Z(u) is a smooth function which yields the corresponding mode number (times tt) 
whenever u equals a hole or magnon root. The mode numbers clearly "label" or "count" 
the solutions of the Bethe equations, and the counting function smoothly interpolates 
between them. 

To write down the one-loop non-linear integral equation, we recall [22] that for an 
arbitrary function f(u), which is analytic within a strip around the real axis, the following 
identity holds 

e/w+e ^ = - r Pf( u ) z ( u )+ r imiog [i + ( -i)v^+*°)] . 

(2.16) 
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Applying this identity to Z{u) and adapting the steps of [32J to the present case, we 



r(l/2 + iu) A., V\-i{u-^] 



Z{u) = iL\og 9 _. t 



r(i/2-i«) ^ r (<(«-«£> 



+ hm f ^z -flog r H ~ i m i og [i + (-1)^*^0)] , (2 . 18) 

a^oo y_ a 7T GSM Y\l [U — V )) L J 

The identity (I2.16P may also be used to express all conserved charges in terms of the 
counting function. The first charge (the momentum), however, needs to be regularized 



P = lim I - 

a— >oo 



r ^ p,{u) z{u) ~ tt p{u ^ + 1- ^ p,{u) imi ° g [i + { ~ i)s eiz{u+io)] 

(2.19) 

In the above formula p(u) denotes the momentum of a single particle 

p(u =-bg 4-. 2.20 

z u — 1/2 

Due to antisymmetry of Z(w) and one easily finds 

P = 0. (2.21) 
Similarly, the one-loop anomalous dimension 71, see ( 12.21) . may be rewritten as 



71 = 4 7E L + 2 ^ {^(1/2 + z u^) +V(l/2-i u^) } 

+ 2 /3 4 (^R^tB^) t 1 + (-1)^-^°)] , (2.22) 

where 7e is Euler's constant. 

Note that the NLIE (12.181) in conjunction with the Bethe equations for the hole roots 
(I2.14p is fully equivalent, for the ground-state, to the algebraic Bethe equations (12.11) for 
arbitrary finite values of M and L. (The generalization to the case of excited states is 
fairly straightforward but will not be discussed in this paper.) Likewise, the expressions 
for the one- loop anomalous dimension 71 given in (12.21) and (I2.22p are equivalent. 



7 Due to superficial divergencies one needs to apply (|2.16[) to Z"(u) and then to integrate twice the 
resulting equation. The constants of integration are fixed by antisymmetry of Z{u) and the condition 

lim Z'{u) = 0. (2.17) 



9 



2.3 Magnon Density 

If the number of magnon roots M gets large we may expect, for the groundstate, that they 
form a dense distribution on the union of two intervals [—b, —a] and [a, b] on the real axis. 
This allows us to introduce a distribution density p m (u), see section 3.2. of [6 J for further 
details. It then follows from (TjLHD and (EIED that 



^Z{u) = 2 vr p m (u) + 2 vr 8{u) + O , with 2 J du p m {u) = 1 , (2.23) 



= 7 7T 0_ 4- 2 7T f Mil) 4- H ( — \ With 2 / 

M dw 

where the ^-function stems from the gap in the center of the magnon mode numbers 
Using this relation one can rewrite ( 12.121) as 

2« Pm{ u) +2^i(«) - Yi7F+\ - 2 (f-° dv+ l dv ) T^WTi = °- (2 ' 24) 

where u G [—6, —a] U [a, 6]. If there is a gap 2a > we therefore may drop the term 
involving the 5-function in ( 12.24ft . In principle, if interpreted appropriately, this equation 
should hold for large M and arbitrary, small or large, L. 

If in addition L stays finite (but arbitrary) we can apply the scaling procedure u = u/M 
of [B], as in this case the gap 2a closes (a — > 0), and in addition b — > M/2. Then the 
non-singular integral equation (I2.24p turns into the singular integral equation, and with 
Po(u) = M p m (u) we find 



Po(u' 



-47r5(n)-2 4 du' / uv \ = 0. (2.25) 
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The solution is the (singular) density 



. , 1 , 1 + V1-4M 2 



- log , = - arctanh y/l -4m 2 ) , (2.26) 

7T 6 1- VI -4M 2 7T V / V 7 

first derived in [37]. It should be considered as a distribution (in the mathematical sense) 
rather than as a regular function. The reason is that the expression for the one-loop 
anomalous dimension in (12.21) formally turns into Am f du p m (u)8(u) = 47rp m (0) = oo. 
However, a more careful analysis [6j of the multiplication of the distributions p m (u) and 
S(u) leads to 

7l = 8 log M + O{M ). (2.27) 

If instead L —>■ oo along with M — > oo such that (3 = M/L is kept finite, the gap 2a 
does not close. We may then drop the 5-function term in (12.241) and obtain after rescaling 
u = u/M with a = a/M, b = b/M 

•—a fb \ i 

= 0, (2.28) 
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which is essentially (up to a rescaling of u by (5) the derivative d/du of the singular two- 
cut integral equation first derived in [TU]. The original equation is easily reconstructed 
by integrating both sides of (12.281) w.r.t. u, with a constant of integration of 2tt//3 on 
the right hand side. The explicit solution for the density p m {v) along with a, b was also 
given in [16J. When (3 — > oo the gap 2a disappears and the limiting distribution (12. 26f) 
is recovered. However, this procedure does not reproduce the correct behavior of the 
anomalous dimension of the previous large M limit at fixed L, i.e. (12.271) ; instead, one 
finds 

7i = | log 2 f + ..•• (2.29) 

See also the discussion in [6]. We notice that large M analysis is quite subtle if the gap 2a 
is very small but non-vanishing. 

In fact, there is a very interesting perturbation on the scaling behavior (I2.27P first 
noticed in [3J. Let us understand this effect by a more refined analysis of (12.231) . (12.241) . 
It is convenient to split the density p m (u) into the singular, leading piece po(u) and a 
fluctuation correction a(u): p m (u) = po(u) + a(u) where po(u) = 1/M p (u/M), see (12.261) . 
The trick is to now add 

2/ dv , P °^] = 4 b f ^ (arctanftt + a) - arctanfn - a)) + O(M ) (2.30) 

J_ a (u — V) 1 + 1 71 M 

to (I2.24H . We then see that <j(u) scales as log M/M and we should therefore define, in 
analogy with |6J, a fluctuation density a{u) through 

8 logM 

Pm{u) = po{u) — — a(u) . (2.31) 



It satisfies 



1 j 1 

2 7T o(u) — — (arctan(-u + a) — arctan(w — a)) + - — -r (2.32) 

2n 8 u 2 + 7 



00 



2 I / r/r- / dr 1 ; a( - V ] = . 



(w-t>) 2 + l 

This integral equation fully determines the fluctuation density a{u), as the edge parameter 
a may be determined from the normalization condition (j_^ a J r / a °°j dup{u) = 1, which 
implies 

j = 8 / aW(w) (2.33) 

The one-loop anomalous dimension is then given from (12 .2j) by 

71W _ G "i .,. ,, ( f " , r , , 



logM 7T 



arctan2a- 16 ^ rfw + jT dnj . (2.34) 
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2.4 Fourier Space Equation 

It is very instructive to change from it-space to Fourier space. After rewriting ( 12. 32ft as 

1 j 1 
a(u) = — - (arctanfu + a) — arctanfu — a)) T (2.35) 

47T 2 16 71 11 -\- -T 



f oo 

+ ' 



dv cr(v) f a dv cr(v) 



IT 1 + [U — V) 2 J_ a TC 1 + [U — V 



and Fourier transformin 



/oo 
due~ itu a(u) (2.36) 
-oo 



one obtains 



&{t) = \K h (t, 0; a) - -L - 4 ^ df ^ h (t, if- a) &{lf)j , (2.37) 

where the kernel is given by 

K h (t,lf;a) = ^t J du cos(tw) cos(t'u), (2.38) 

which leads to the expression ( II. 151) in the introduction. Likewise, Fourier-transforming the 
normalization condition (12.331) yields the relation (11.161) between the physical parameter j 
and the fluctuation density a(t) in Fourier space stated already in the introduction. The 
one-loop anomalous dimension is then given by 



logM 



2 

1 arctan2a (2.39) 

7T 



7i(i) 

- = o ± - 

poo / poo \ 

4 J dt(a{t)-At J dt' k h (t,tf;a) a(tf)J 



2.5 Hole Density 

The Bethe roots corresponding to the small holes lie inside some interval [— c, c]. In the 
"thermodynamic" limit L — > oo, where the number of small holes tends to infinity, their 
one- loop root distribution density Ph(u) is related to the counting function through 



\^Z{u) =27rp h (w) + Qrj , with J du Ph {u) = l, (2.40) 



8 We have included a factor of e 2 into this definition for convenience. For all other Fourier transformed 
quantities in this paper, in particular all kernels K, we do not include such a factor. 
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as one easily derives from (I2.1(jp and (j2.14p . Using (12.181) . we may then derive a nonlinear 
integral equation for the distribution of holes 

p h (u) = j (u - + j>(-i (u - + ^(i (« + + V>H ( u + 



r-c 

+ 



jf ^ (it - „)) + Vi(-< (« - «))) 0h(«) , (2.41) 



where the term -| j-I(u) denotes the derivative of the last line in (12.181) . The terms on the 
r.h.s. of the first line of (12.411) are the contributions of the two large holes with rapidities 
u^\u^ = — «k j cf (12.111) . where we have also implicitly assumed L <C M. Then the two 
rapid holes behave aa^| u^' 2 ^ ~ ±M/ ^/2, while the term 1/L -^X{u) in (12.411) yields merely 
an additive 2 log 2, see appendix |X] for a discussion of this point. The four terms on the 
r.h.s. of the first line of (12.411) thus behave like 4 log Mj a/2. Using (11.71) we thus derive a 
linear integral equation 

2 1 f c dv 

PM = — {^+iu)+^(i-iu))+ J — (i>(i(u-v))+^(-i(u-v)))p h (v). (2.42) 

One then finds the generalized one-loop scaling function, cf. (II. 8p . from (12.221) 

~' lU) °» + 2j I dup h {u)(^{\ + iu)+i){\-iu)-2^{l)). (2.43) 



logM 



In order to easily generate the series expansion of (12.431) in powers of j , defined in (11.71) , 
it is useful to rescale u and define 

u 

u = - and p h (u) = j cp h (u). (2.44) 
Defining the non-singular kernel 

K(u, v) = y4>{i c(u — v)) + ip(—ic(u — v)) — + icu) — — icu)j , (2.45) 

the integral equation (I2.42p becomes 



2 



p h (u) = -c+ dvK(u,v)p h (v) . (2.46) 
n J-i 

It is of Fredholm-type and may be immediately expanded in the small parameter c and 
iteratively solved as a power series in c. The relation to the parameter j is then determined 
through the normalization condition in (12.401) which becomes 

j= [ dvp h (u). (2.47) 



9 Extensive numerical studies indicate that for the ground state at large M all charges qi in (|2.5p are 
small except qi- Then one finds from t(u) = and §F]\ u£' 2) ~ ±M/\/2. See also 0. 
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This yields j as a series in c. The generalized one-loop scaling function (I2.43P becomes 

j^£: = 8 + 2 J dup h {u){ij{\ + icu)+^{\-icu)-2ij{l)). (2.48) 

This yields the one-loop scaling function as a series in c. Inverting the series (12.471) and 
substituting into the expansion of (12.481) gives the desired series of the scaling function in 
terms of j. It starts out as 

-^4t = 8-8j log2 + — j 3 7r 2 C(3)--j 4 7r 2 log2C(3) 
logM J 6 12 J sv ' Q J 6 sv ; 

+ 2f (^n 2 log 2 2 ((3)--^^ ((5)^ +0(f). (2.49) 
Note that by analytic continuation the density of the holes is related to u{u) via 



2 

j p h (u) = 8a(u) ue(-c,c) (2.50) 

7T 



which may be rewritten as 



2 8 

j Ph(u) = / dt a(t) e* costu . (2.51) 



7T 7T 



The preceding derivation proceeds from the counting function, cf (12.401) . We will closely 
follow this procedure in the next chapter [31 where we will treat the higher-loop case. It 
should be noted, however, that our solution (12.461) . ( !2~471) . (12.481) may also be immediately 
recovered by Fourier analyzing the results of the previous section 12.41 The reader should 
multiply (12.371) with ea costu, integrate in t over the positive real semi-axis and use the 
integral representation of the kernel (I2.38p . Subsequently rewriting j in terms of a(t) with 
the help of (jl.ISp and finally using the relation (I2.5ip it is straightforward to derive (12.421) . 
We thus conclude that 

a = c. (2.52) 

This equation tells us that the gap [—a, a] in the distribution of magnon roots is densely 
filled by the (small) hole roots. 



3 All-Loop Theory 

3.1 The Asymptotic Non-Linear Integral Equation (NLIE) 

Let us now extend the one-loop results of the last chapter to the higher loop case. We will 
use the asymptotic Bethe ansatz for AdS/CFT, based on the S-matrix approach [38]. In 
the sl(2) subsector the asymptotic all-loop Bethe equations j5[ [7] read 

M . I 1 — 



^ = TT ^~^~ Z g£ e 2i ^\ (3.1) 
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We define the all-loop asymptotic counting function as 

/ • In \ M 

f s. . rl x{i/2 + u) . . y^i i + u-u k 
Z{u) — %L log ^rj^ + , £ log - (^) 

M i _i 2! 

x(i/2+u)x(i/2-u k ) 



_ 2l J2\og ^V2+uMi/2- Ufc) + ^ ^ (3 2) 

fc=l x(i/2-u)x(i/2+u k ) k=l 

As in the one-loop case, one finds the corresponding non-linear integral equation 

— 0'(u - v, 1) Imlog [1 + (-1) 5 e iZ ^°)] 

-co " 



L 




j'=i 




POO 




J —00 


2^ 


L 





"CO 

.(?) 



-00 

„2 



I I , , • ^ 1 x(i/2+u)x(i/2— v) a 1 v 

-/ ~ I 2 ^ lQg i I ^ 

' s(i/2-«)a!(i/2+«) 

x(i/2+u)x(i/2-u^) . 



£ 2 i log - ggggggiJ _ 9{U) u t 

\ x(i/2— up 

dv I '. d , 1 + 



2 .^_ bg _ g( i/2+uMi/2-,) _ g ^ >7;) j imlog [l + (-l)-5 e ^(^o)] 



x(i/2— u)xti/2+v) 



(3.3) 



The counting function defined in (I3.2p satisfies a similar relation to f)2.23p . but with the 
all-loop density on the r.h.s. 

3.2 The NLIE in Fourier Space 

In Fourier t-space equation ( 13. 3ft becomes 



o„ r„| £ 2tt cos o 

w it(e*-l) ^ it(e*-l) e* - 1 

+ 8<? 2 ^— y ^ dt' e~2 K(2gt,2gt') ^ C(t') 

L 

71 x - 
+ -> • 

i=i 

-4^ 2 ^— / dt'e--t'K(2gt } 2gt')Z(t'), (3.4) 
e — 1 Jo 



^ ^ cos ( t' v!jp 

1 3=1 
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where C(t) denotes the Fourier transform of the "Im log" term. Note that the Z(t) has a 
first order pole at t = 0. This in accordance with (13.21) . since the Fourier transform of this 
expression must be understood in the principal value sense. Note that we have not made 
any approximations. Therefore H3.4I) is still fully equivalent to the orginal set of discrete 
asymptotic equations ( 13. ip . 



3.3 Large Parameter Integrals 

Let us now investigate the effects of taking the large M limit with L <C M. It will be 
important to understand the large M expansion of integrals of the form 

POO 

f(M)= dxh{x) sm{u{M)x) , (3.5) 



where h(x) is a smooth integrable function on [0, oo) and u(M) — > oo when M — > oo. 
We first note that because of the relation to the Fourier transform (Plancherel's theorem) 
limjv/^oo f(M) = 0. Since f(M) is meromorphic and vanishes at infinity we have 

oo 

To find Co it is sufficient to note that 

c = lim u(M) f(M) = lim / dx h(x) ( — ^ cos (u(M) x) ) = h(0), (3.7) 
m^oo m->oc J \ dx J 

since the integral after a partial integration vanishes again. By subsequent integrations by 
part one finds that 

c n = lim u(M) n+1 ( f(M) - V — — ^p— — ] = (-1)5 h {n \0) for even n . (3.8) 

M— >oo y ^ U\M) L+:) J 

The odd c n coefficients vanish, as follows from (I3.5p . 



3.4 The Leading Order Equation 

To derive from (13.41) an equation reproducing the leading contribution to the scaling func- 
tion in the limit where M — > oo and L is kept fixed, it is sufficient to observe, based on 
the results of the previous subsection, that upon iterating ( 13.41) only terms of the form 



(1,2) 
h 



27re | 2vrcos(tw 
it (e* - 1) it (e* - 1) 



(3.9) 
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where u^' 2 ^ ~ ±y|g 2 ~ represent the universal holes, will give the leading (loga- 
rithmic) contribution. This is because we have 

u^~0 j = 3,...,L (3.10) 

at leading order, and the terms involving C{t) do not contribute at this order, see appendix 
lAl Thus the leading all-loop equation reads 



4vret 4 7 rcos(t4 1) 
Z t) - v 



it(e*-l) it(e*-l) 

-4^ 2 ^- / dffe~TlfK(2gt,2gif)Z(1f) (3.11) 
e* - 1 Jo 

Upon subtracting the one-loop part of this equation 

Z(t) = Z (t) + 5Z BES (t) (3.12) 
and identifying 8Z(t) with the fluctuation density 

5Z BE3 (t) = Wmg 2 e l * log(M) (3.13) 

one rederives the equation of [?]: 

<W*) = (k(2gt,0)-4g 2 J dt' K{2gt,2gt') <7 BES (t')) . (3.14) 

3.5 Subleading Corrections to the Twist Operator Dimensions 

The large M expansion of the anomalous dimensions of twist operators is expected to have 
the following form 



7 = f{9) logM + f sl (g, L) + [-^y ) , (3-15) 

where f s \(g,L) denotes the subleading effects of O(M ). These are easily obtained from 
(13. 4p . and we may compute f s \(g, L) to arbitrary order of perturbation theory: 

f sl (g,L) = ( 1 -(L-2) log 2) f(g) - 8 (7 - 2 L) C(3) g 4 
+ § i^A C(3) + (62 - 21 L) C(5)) g 6 



-— ((13-3L)7r 4 C(3) + 5(32-llL)7r 2 C(5) + 75(127-46L)C(7)) g 8 
±... (3.16) 

Notice that the "universality", i.e. L-independence of the scaling function f(g) is lost when 
one computes these O(M ) terms. They contain L-independent and terms linear in L. 
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4 The Generalized Scaling Function 



4.1 Derivation 

Let us now finally treat the novel scaling limit (11. 7ft . i.e. we consider the limit L, M — > oo 
with j = L/logM kept fixed. In this limit, in contradistinction to section [331 a ^ so the 
L — 2 remaining holes contribute. Although individual hole terms separately do not develop 
logarithmic behavior in M, their collective contribution will be proportional to L = j log M. 
Furthermore, in this limit all terms involving C(t) can be dropped, see appendix [XJ Thus 
( 13. 41) for the counting function Z(t) in Fourier space linearizes in this limit to the form 



Z(t) 



2nLe 



it(e 



L 2 7T COS ( t U 



,00 



zt(e* - 1' 



L-2 



+ 8tt g 1 



z (e 

t 

e-2 



>' - 1) <^ 



dt' e" K(2gt, 2gt') cos t! u( 



,(i) 



1 



dt' t' K(2gt,2gt') Z[t') 



(4.1) 



Note that in above formula only quantum corrections to u J h for j = 3, . . . , L need to be 
taken into account, since the corrections to the universal holes are, upon the iteration, 
subleading. In similarity to section 13.41 we strip off the one-loop part by defining 



Z(t) = Z (t) + SZ(t). 
We relate SZ(t) to the fluctuation density a(t) through 

bZit) = 16me^ logM. 



(4.2) 



(4.3) 



and derive to the desired order 
t 



a{t) 



e* - 1 



j Jo(2gt) 



+ 



~*/ 2 cos(t 



81ogM ^ 



2 i n 

9 1 V- 

2 logM ^ 

6 i=3 



dt' K(2gt, 2gt') e"*' /2 cos(t' u£ 



(4.4) 



-4^ 



dt' K(2gt : 2gt')a(t') 



The corresponding anomalous dimension can be easily shown to be given by 

L 



j = 8g logM 



logM ^ 



o 



2gt 



cos(tu 



h / 



(4.5) 



Ji(2^) 
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The distribution of the small holes is found from 

Z(u{) = n(2n{ + 6-l), (4.6) 

which in Fourier space reads 

i 

Plugging (14. 2 p into (14. 7\\ and observing that (see section 13 . 3f) 

t 

e2 — cos ty 



oo 

sin (tu{) Z(t) = n(2n{ + 5- 1). (4.7) 



POO 

F'(x,y) = / dt costx 



= ~ ^ (i (x - ?/)) + ip (-i {x -y))+ip (i (x + y)) + ip {-% (x + y)) 

-2^Q-ix) -2vQ + ix^ Y (4.8) 

one easily derives from ( 14.71) 

2 7in { h k) =4F(u { h k \u^) - 16 logM f°° dt^-e t/2 sm(tu^). (4.9) 

Jo t 

Introducing the density of holes Ph_{u) it follows from (14. 9p that 

Q POO 

, .■F'fitVV- / dta(t)e t/2 cos(*u). (4.10) 
vrlogM v h ; vr y w v ; v ; 

Note that - A F'(u, ) is at large values of M essentially the Korchemsky density po(u), 
i.e. (12.261) after scaling back u = M u, po(u) = 1/M p (u), up to small corrections at the 
boundaries of the distribution of the roots. Since the small holes occupy a finite interval 
(—a, a) one can safely take the large M limiS 

F'(u,u^) = log M + O(M ) ue(-a,a). (4.11) 

After replacing the sum in (14. 4p by an integral and using the above density we find 
t 



( \ 2 
3Ph{u) 



a{t) 



-j-J (2gt) + K h (t,0;a)-4 [ dt' K h (t,t';a)a(t') 
. 8t Jo 

POO 

+g 2 K{2gt, 0)-4g 2 / dt' K{2gt, 2gt!) a{t') 
Jo 

/oo / poo 

dt' t' K(2gt, 2gt') ik h (t', 0; a) - 4 J dt" K h (t', t") a{t"\ 



(4.12) 



10 The magnon density is related to a(t) by a similar formula 

2 1., m, 8 1ogM f 00 , , , 
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where Kt(t, t';a) is the one-loop kernel given in (j!.15p . The endpoints can be obtained 
from the normalization condition 

pa 

dup(u) = l (4.13) 

which implies ( 11.161) . Inserting (11.161) into (14.121) we find the final integral equation (ll.9p 
anounced in the introduction. Likewise, the anomalous dimension (14.51) may be reexpressed 
and simplified as 



7 = 8 g 2 log M 



dt^M^ h (t ,0; ;a) 
2gt v ' ' ' 



8 f°° dt Jl ^ 2g ^ (a{t) -At r dt'K h (t,t';a)a(t')) 

J0 2(yft \ JO J 

16 logM (a(0) + 1- J . (4.14) 



16, 

This concludes our derivation of the equations determining the generalized scaling function 
f(g, j) in ( 11.81) . Let us now apply them to obtain the first few terms in the double expansion 
of this function in powers of g and j. 

4.2 Weak Coupling Expansion 

The equation (14. 12ft is solved iteratively with relative ease in a double-perturbative series 
in g and the gap parameter a. As in the one-loop case in section [2] one then inverts (11. 16ft 
to obtain a(j) as a power series in j. This then yields the fluctuation density a(t) as a 
series in g and j. It starts out as 

a(t) =g 2 a BEa (t) 

1 1 g 2 t (t - 4\og2) 



+ J 



8 e*/ 2 + e* 8 e* - 1 



+ # 4 — (-3t 3 - 4vr 2 t + 16vr 2 log2 + 24t 2 log2 + 96C(3)) + ... J 
— 1 96 J 



+ f x 



.3/ 7T 2 2 t , //on #71" t (14£(3) — it te / ) 

+ 7 I tVcsch t/2 + - v sv ^ i 

V !536 w y 384 e* - 1 

+ ^ - 1 (37r 4 t - 7r 4 te-*/ 2 
2304 e* - 1 V 

+ 140tt 2 C(3) - 42C(3)t 2 - 2232C(5)) + . . . J 
+ ... (4.15) 
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The generalized scaling function at weak coupling is simply given via (11 . IT)) by evaluating 
the fluctuation density at t = 0. Let us define an infinite set of functions {f n (g)} as 

oo 

f(9j) = f(g) + Y,fn(g)j n - (4.16) 

71=1 

The first one fi(g) is 

h (g) = -8 g 2 log 2 + g 4 Q vr 2 log 2 + 16 C(3)) - g< vr 4 log 2 + ~ rr 2 C(3) + 168 C(5)) 
+ g 8 (j^ log 2 + | tt 4 C(3) + 64 log 2 C(3) 2 + f tt 2 C(5) + 1840 C(7)) + . . . 

(4.17) 

Note that /i(<?) is special as it can be obtained from ( 13. 161) by keeping only terms propor- 
tional to L. To this order the hole momenta are set to zero. Only at orders higher than 
linear in j one needs to take into account the "dynamics" of the holes. We then find for 

h(g),---U(g) 

fi(g) = -fig) log 2 + 16 g 4 C(3) - g 6 Q vr 2 C(3) + 168 C(5)) 
+ g 8 tt 4 C(3) + f rr 2 C(5) + 1840 C(7)) + . . . 

(4.18) 



/ 2 (ff) = 



(4.19) 



M = y 2 9 2 * 2 C(3) + g 4 (§ vr 4 C (3) - I vr 2 C (5) 



+ .. 
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540 
7 



^ C(3 ) + 1^ C(3) 3 + - tt 6 C(5) + — vr 4 C(7) 
108 sv ; 3 sv ; 15 sv ; 6 SV ; 



17885 



tt 2 C(9) 



(4.20) 
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7 / 77 4Q N 

U{g) = ~q9 2 * 2 log 2 C(3) + g 4 (-- vr 4 log 2 C(3) + - tt 2 C(3) 2 + 31 vr 2 log 2 C(5) 

+ 96 (- log 2 c(3) + ft " 4 c(3)2 + x " 4 log 2 c(5) 

fi51 

-— vr 2 C(3)C(5)- 635 7r 2 log2C(7) 
+ 9 8 tt 8 log 2 C(3) + g tt 6 C(3) 2 - 252 vr 2 log 2 C(3) 3 + ^ tt 6 log 2 C(5) 
' :,S)1 1 vr 4 C(3) C(5) + 2883 tt 2 C(5) 2 - ^ tt 4 log 2 C(7) 



U 

17780 2 > , , >/fT , 35770 2 , 
+ ^— vr 2 C(3) C(7) + —3- tt log 2 C(9) 

+ .... (4.21) 

At fixed j, we observe a constant degree of transcendentality [39] of all terms contributing 
to a given order of perturbation theory in the coupling g. Interestingly, the converse is not 
true, as may already be seen from the one- loop result f)2.49p . 

As was announced earlier the function /2G7) is identically zero, indicating that all 
terms of order j 2 in the j-expansion of f(g,j) are absent to all orders in the coupling 
constant g. This is easily proven directly from our equations. Some potentially related 
very interesting observations at strong coupling were made in [23]. Roiban and Tseytlin 
found some intriguing evidence that terms of the form j 2 \og k j might upon resummation 
indeed result in a vanishing j 2 contribution, cf also the discussion in the introduction. 
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A The Non-Linear Term 

In this appendix we will discuss the integrals involving the non-linear term. For simplicity 
we will confine ourselves to the one-loop case, where it is sufficient to consider 

J(«)= lim ^ - ^log T( \^ U ~ f Imlog [1 + (-l)' e «*(«*o)] (Al) 

We first note that the function 

C{u) = Imlog [1 + (A.2) 
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is smooth apart from a finite numbers of points, namely when u is equal to the magnon or 
the hole rapidity. A closer inspection reveals that 



C{ui — e) = 7r C{ui + e) = — n (A. 3) 

where Ui denotes either a hole or a magnon rapidity We will assume that the small holes 
and the magnons are densely distributed along the real axis, as this is the case for the 
limits discussed in this paper. It is easy to convince oneself that the integral (1 A. II) gets 
the dominant contribution from (—a, — 4^-) U (4f , ol). Because the small roots and magnons 
are, at large values of M, densely and symmetrically distributed on (— 4^) this part of 
the integral contributes starting at O f-rfs) only. Assuming v £ (—a, — 4f) U(4f- , a) we may 
expand the integrand in a power series in u. Because of the antisymmetry of the counting 
function only odd powers of u survive the integration. Thus we may write 

^ lQg r r^-",)) ) = 1 (M-iv)-Miv))u - 1 (u-iv ) 

+ 0(u 5 ) + even terms in v . (A. 4) 

On the other hand from the definition of the counting function we have 

„, . L + 2M ( 1 \ M .. . 

£( v ) = n + v > —F ■ ( A - 5 ) 



2v \v 3 J y/2 

Plugging (1A.4I) and (lA.5j) into (1A.1I) we find 



it 3 



X(u)=tu + 0[^). (A.6) 

To fix the constant ^ it is necessary to extend the expansion in (IA.50 to the whole interval 
v £ (^,00). However there is a much simpler method. Since the above discussion is not 
sensitive to the value of L we may set L = 2. Then we may compute the corresponding 
anomalous dimension plugging 1 12.181) together with ( 1A.1I) into 112. 16f) . Comparison with 
the exact one- loop result 71 = 8 Si(M) fixes £ to be 

£ = 2 log 2 (A.7) 



Numerically we have checked that the expansion 1IA.6I) breaks only around the small neigh- 
borhood of ±4p This suggests that the radius of convergence of 1IA.6I) lies closely to the 
edge of the magnon distribution. 
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